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The complete 2-loop QED contributions to the muon lifetime have been calculated 
analytically in the Fermi theory. The exact result for the effects of virtual and real 
photons, virtual electrons, muons and hadrons as well as e+e — pair creation is 
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where To is the tree-level width. This eliminates the theoretical error in the ex- 
tracted value of the Fermi coupling constant, Gp , which was previously the source 
of the dominant uncertainty. The new value is 

Gp = (1.16637 ± 0.00001) x 10~ 5 GeV -2 

with the error being entirely experimental. Several experiments are planned for the 
next generation of muon lifetime measurements and these can proceed unhindered 
by theoretical uncertainties. 



1 Introduction 

The three fundamental input parameters that enter into all calculations of 
electroweak physics are the electromagnetic coupling constant, a, the Fermi 
coupling constant, Gp, and the mass of the Z° boson, AL^ Their current best 
values, along with their absolute and relative errors areEm 

a = 1/(137.0359895 ±0.0000061) (0.045 ppm) 

G F = (1.16639 ± 0.00002) x 10" 5 GeV" 2 (17 ppm) 

M z = 91.1867 ± 0.0021 GeV (23 ppm) 

The first of these, a, comes with the liability the 'hadronic uncertainty' arising 
because, when used for the analysis of data near the Z° resonance, it must be 
run up from a scale q 2 = to Af§ crossing, on the way, the hadronic resonance 
region. 

In the mid-80's, just before the turn on of LEP, a CERN report concluded 
that the error on Mz would be ±50 MeV or 550 ppm and that "A factor of 2-3 



improvement can be reached with a determined effort" a. It was thus generally 
believed that the error on Mz represented the limiting factor in the precision 
with which theoretical predictions could be made. The situation has changed 
adiabatically, however, and the relative error on Mz now approaches that of 



Since LEP, as a machine, is not a radically new design, the lesson that we 
should take is that it is extremely difficult to predict the accuracy with which 
physical quantities will be measured, even in the relatively short term, and 
that one should constantly strive to reduce such errors to the minimum level 
consistent with the available technology. The possibility of precision physics 
at a muon collides serves to emphasize this point. 

With this in mind, and given the great cost and effort that was expended 
in reducing the error on Mz to its current value, it is reasonable to look again 
at Gf and see what is required to reduce its error to a level where it can never 
become an obstacle limiting the accuracy with which theoretical predictions 
can be made. 

Gp is extracted from the measured value of the muon lifetime, r M = 
(2.19703 ± 0.00004) ^s El and on the experimental side this is currently the 
source of the dominant error. New experiments are planned at the Brookhaven 
National Laboratory, the Paul Scherrer Institute and the Rutherford- Appleton 
Laboratory and it is likely that the uncertainty on Gf from this source will be 
reduced to somewhere in the range 0.5-1 ppm. 

Most of the work reported here appears in ref.s |||(|. 

2 The Fermi Coupling Constant 

As given above the current relative error on the Fermi constant is SGf/Gf — 
1.7 x 10~ 5 . Of this 0.9 x 10~ 5 is experimental and 1.5 x 10~ 5 is theoretical 
being an estimate of unknown 2-loop QED corrections. 

is related to the measured muon lifetime, r, by the formula 
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as calculated using the Fermi theory in which the weak interactions are de- 
scribed by a contact interaction. Aq encapsulates the higher order QED cor- 



rections and may written as a perturbation series in a r = e^/(4n), the renor- 
malized electromagnetic coupling constant. Thus 

oc 
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in which the index i gives the powex. of, a r that appears in Aq^ . Note that 
Eq. differs from the usual formula!! in S3ys that begin to become important 
at the part-per-million level. It is knowrQEa that 

A<? (0) = -8.t - 12.T 2 In x + 8x 3 - x 4 (4) 

Aq{1) = (v) (f " 3C(2) ) + °^ xlnx ) ( 5 ) 

where x — m^/m 2 and C, is the Riemann zeta function with £(2) = 7r 2 /6. 
That the Aq^ remain finite in the limit m e — > is a consequence of the 
Kinoshita-Lee-Nauenberg theorem] whose discovery was largely prompted by 
this particular observation. ._. 

Although the Fermi theory is not renormalizable, the Aq^ can be showrciJ 
to be finite for all i. This remarkable feature follows from the fact that the 
V — A interaction is invariant under a Fierz rearrangement that interchanges 
the wavefunctions of the electron and the muon neutrino. Thus Fermi theory 
is equivalent to an effective theory in which the muon and electron occupy 
the same fermion current in the weak interaction lagrangian. After fermion 
mass renormalization is performed the divergences in the vector part of this 
current are independent of the fermion mass and hence cancel in exactly the 
way they would for QED. The lagrangian of Fermi theory is invariant under 
the transformations ip e — > 75"0e and m e — > ~m e . The QED corrections to 
the axial vector of the part can thus be obtained from the those of the vector 
part by changing the sign of the electron mass and hence are finite as well. 
Moreover the two sets of corrections are equal in the limit m e — > and, in 
that case, calculations need only be performed using the vector part of the 
Fermi interaction. This conclusion holds under any regularization prescription 
and avoids the complications associated with the use of 75 in dimensional 
regularizatioiO. 

The foregoing discussion does not apply to the /3-decay of the neutron 
where the Fierz rearrangement generates scalar and pseudoscalar terms that 
bear no resemblance to QED and the radiative corrections are consequently 
not finite. 



3 The 2-loop QED Corrections to the Muon Lifetime 



The complete 2-loop QED corrections to the muon lifetime require the cal- 
culation of matrix element for the processes, fi~~ — > v^e~v ei ji~ — > v^e~v e ^i 1 
fi~ — » f M e~i> e 77 and — > v^e~i> e e + e~ with up to two virtual photons. All 
processes contain infrared (IR) divergences coming from either virtual pho- 
tons, soft bremsstrahlung or both. The cancellation of IR divergences occurs 
between the various processes but this complication may be avoided by ex- 
ploiting cutting relations and calculating the 2-loop corrections as imaginary 
parts of 4-loop diagrams, some of which are shown in Fig.s 1 and 2. In these 
Fcynman diagrams thick lines represent a muon and the thin lines represent 
either the electron or the neutrinos all of which are taken to be massless. Since 
the external muon is on-shell any cut passing through a muon line will vanish 
and the only cuts contributing to the imaginary part are precisely the ones 
that generate the diagrams appearing in the calculation of muon decay. 

Recursion relationslia obtained by integration-by-parts were first applied to 
reduce all dimcnsionally regularized integrals to a small set of relatively simple 
integrals. The well-behaved primitive integrals were then calculated by taking 
the external muon momentum, q, off mass shell to obtain expressions as power 
series in x = — 3 /miLand logarithms of x using well-established large mass 
expansion techniques!^. As the large mass expansion proceeds many terms, 
such as those that are topologically tadpoles, can be immediately discarded 
since they do not give rise to imaginary parts. Since the final result is required 
for x = 1 the complete series must be summed which can now be done in closed 
form in terms of polygamma functions and certain classes of multiple nested 
sumsEi. 

All diagrams were calculated in a general covariant gauge for the photon 
field and exact cancellation in the final result of the dependence on the gauge 
parameter was demonstrated. 



3.1 Photonic Corrections 

Examples of photonic diagrams which when cut give rise to contributions to 
the muon lifetime at 0(a 2 ) are shown in Fig. 2. 

The result obtained for the complete set of photonic diagrams is 
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Figure 1: Examples of diagrams whose cuts give contributions to fi 
^e _ P e 7 or fj,- -> ^e-p e 77. 



where C(3) = 1.20206... and C(4) = vr 4 /90. 
3.2 Electron- Loops and e + e~ Pair Creation 

Diagrams containing an electron loop whose cuts give contributions to muon 
decay are shown in Fig. 2. The result obtained for these diagrams is 
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The value-given in Eq. ([|) is consistent with a numerical study carried out 
by Luke et a/cJ in the context of semi-leptonic decays of heavy quarks. 

In order to obtain a UV finite answer the a diagrams in which the electron 
loop is replaced by the photon 2-point counterterm must be included and 
therefore a decision has to taken as to the renormalization scheme that is to 
be adopted. This will be discussed further in section 4. 



3.3 Hadronic Contributions 

Hadronic effects enter at the 2-loop level through the diagrams shown in 
Fig. 3. The shaded blob represents the hadronic vacuum polarization of the 



photon. The hadronic contribution can be calculated in the usual way using 
dispersion relations but, in contrast to other well-known situations for which 
such effects have been calculated Ej'EZI, the momenta of the external fermions, 
to which the virtual photon is attached, is not fixed. Here the electron par- 
ticipates in the phase-space integration which complicates matters somewhat. 
Overall the shift induced in the inverse lifetime, r M , of the muon is given as a 
convolution integral 

Ar Uad = ^ [°° -R(mlz) AT(z) (10) 
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over the hadronic spectrum, R{q 2 ) = Uhad/opoint, anf l m which P = m 'i/ rn ft = 
1. 61395.. _ The convolution kernel, Ar(z) is obtained exactly as an analytic 
function B. When the integral is performed using actual hadronic data the 
result is 

-^J (0.042 ± 0.002) (11) 

which includes a rather conservative estimate of the hadronic uncertainty. Still 
the latter amounts to only 2 parts in 10 8 and so is well under control. 

The integral ( JlOj ) can be used to obtain an expression for the contribution 
from diagrams where the hadronic vacuum polarization has been replaced by 
muon loop by setting 

R(mlz) = (l + f ) \f^\- (12) 

which gives 

^-r.£)-(^->->) (13) 

= -r ( — J 0.0364333. (14) 

The result agrees with that obtained by perturbative methods. The effect of 
tau loops can be obtained in a similar way and, as expected on the basis of 
the decoupling theorem, is negligibly small. 

4 The Renormalized Electromagnetic Coupling Constant, a r 

The use of dispersion relations to calculate the hadronic and muon loop con- 
tributions in the previous section naturally invokes a subtraction of the pho- 
ton vacuum polarization at q 2 = and is therefore equivalent to the on-shell 



renormalization scheme. In cases where there are two or more widely sepa- 
rated scales, such as m e and m^, use of the MS renormalization scheme is 
indicated since it automatically incorporates the large logarithms that arise 
into the value of the renormalized coupling constant, a r , at tree level. 

It is therefore appropriate here to adopt the MS renormalization scheme. 
The hadronic contributions of section 3.3 that were obtained via dispersion 
relations must be corrected to convert them from the on-shell to MS renor- 
malization scheme. As it turns out the contribution from muon loops is the 
same in both schemes when the 't_Hooft mass is taken set to fj, = as is 
appropriate here. It can be showiES that the MS renormalization scheme is 
implemented in a a consistent manner by using the results of section 3.3 as 
they are given and setting 




where the logarithm of 0(a 3 ) was first calculated by Jost and Luttingc £l The 
substitution ( fl5l ) correctly resums logarithms of the form a" In™ -1 {m 2 /m^) for 
all n > and incorporates those of a 3 ln(m^/mg). Upon evaluation 

a e (m /1 ) = 1/135.90 = 0.0073582. 

5 Conclusions 

It has been over 40 years since the 1-loop QED corrections to the muon lifetime 
were calculated. The 2-loop contributions have had to await the development 
of new theoretical techniques, as well substantial increases in computer speed 
and storage capacity, but are now available. 

The complete 2-loop QED contribution to the muon lifetime in the Fermi 
model may be encapsulated in the quantity Aq( 2 \ as defined in Eq.s(Q) and 
(||) , and is given by 

= /oeWy / 156815 _ 1036 _ 895 67 
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+ 53C(2) In 2 - (0.042 ± 0.002)^. (16) 
This translates into a new value for the Fermi coupling constant of 



G F = (1.16637 ±0.00001) x 10~ 5 GeV -2 



(9 ppm) 



The error has been halved relative to its previous value and is now entirely 
experimental. 

New measurements of the muon lifetime are planned at the Brookhaven 
National Laboratory, the Paul Scherrer Institute and the Rutherford- Appleton 
Laboratory and it is therefore likely that the uncertainty on Gf from this 
source will be reduced to somewhere in the range 0.5-1 ppm. 

In that case the theoretical error should still be negligible but other issues 
such, as error on the muon mass, m M , and the upper limit on the muon neutrino 
mass, to„ , need to be considered. 

Finally many of the results and techniques employed here can be readily 
taken over and applied to inclusive decays of the 6-quark. 
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